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, ^ Abstract 

Ln ; 

We consider a quantum particle in a periodic structure submitted 
, to a constant external electromotive force. The periodic background 

<^ ' is given by a smooth potential plus singular point interactions and 



has the property that the gaps between its bands are growing with 
the band index. We prove that the spectrum is pure point-i.e. tra- 
jectories of wave packets lie in compact sets in Hilbert space- if the 
Bloch frequency is non-resonant with the frequency of the system and 
satisfies a Diophantine type estimate, or if it is resonant. Furthermore 
it is shown that the KAM method employed in the non-resonant case 
produces uniform bounds on the growth of energy for driven systems. 



1 Introduction 



We study stability of the dynamics of one electron in a Id periodic structure 
with infinitely many open gaps driven by a constant electromotive force. 
To be specific we consider two realizations: the Stark- Wannier problem for a 
periodic background interaction V(x) = V(x + L) defined by the Hamiltonian 



H s = -A + V{x) + Fx on L 
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and an electron on a conducting ring in the plane threaded by a linearly 
increasing magnetic flux line $(t) = Ft whose dynamics is defined by 

H R (t) = (-id x - Ft) 2 + V(x) on L^S 1 ). 

The stability of a general time dependent system was addressed for exam- 
ple in [[H], |K], [13], ^6], [22], [12], ^] . The discussion of stability may be summarized 



in the question whether a wave packet can get delocalized during its time 
evolution. To answer this one may study the time behavior of expectations 
of observables; if the system is periodic in time the spectral properties of 
the Floquet operator-i.e.: the evolution over one period- can provide precise 
information on the stability. To mention one example: if the periodically 
time dependent system is confined and unbounded and the spectrum of its 
Floquet operator is absolutely continuous then the energy expectation grows 
in time for any initial state, so the system is unstable. 

The special case considered here was intensively studied since Wannier 



conjectured existence of ladders of eigenvalues; see [19], for background 
on this story. In [0] it was proven that for smooth background potential 
V G C 2 (R) -in fact C 1+e - the spectrum of Hg is absolutely continuous which 
leads to unbounded growth of the energy for H R (t), see M. On the other 



hand in ||, |L8| a comb of 5' point interactions was considered. It was shown 
that this model is physically important, in particular it describes idealized 
geometric scatterers. It was proven that the spectrum has no absolutely con- 
tinuous component leaving the possibility of eigenstates, singular continuous 
spectrum and unbounded energy growth. Furthermore a conjecture on the 



essential spectrum was made. See also [[23[ for the geometric scatterer aspect 



and |24] for a second proof of absence of absolute continuity. 

It was argued by Ao [|TJ that the spectral nature depends on the gap 
structure of the periodic background. He conjectured that for gap behavior 
AE n = 0(l/n a ) one has point spectrum for a < at least for "non- resonant" 
F and continuous spectrum for a > 0. For a = a phase transition from pure 
point to continuous spectrum with growing F is expected -see also 0, ||; 
furthermore the spectral nature seems to depend also on number theoretical 
properties of the driving frequency FL. This critical case corresponds to 
the driven Kronig-Penney model; another realization of constant gaps is the 
-explicitly solvable- forced harmonic oscillator considered by |17], IT[ [10 |. 



Our contribution, here, is to show stability for V a comb of 6' interactions 
plus a smooth bounded background. We prove that the spectrum of H s 
is pure point; H R is periodic up to a gauge transformation, the Floquet 
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Hamiltonian of the transformed problem is unitarily equivalent to H$, its 
spectrum is also pure point. We are able to prove this in two different settings. 
Firstly for a large subset of frequencies FL which do not resonate with the 
frequencies of the background we use a KAM algorithm in order to treat the 
8' interaction as perturbation of the decoupled problem where 8' is replaced 
by a Neumann boundary condition. This algorithm needs as input a matrix 
which has sufficient off diagonal decay. Because of the singularity of the 
interaction it is not evident that the Floquet Hamiltonian considered here 
has such a matrix representation. A detailed spectral analysis is necessary 
to show that this is the case. Technically the basic observation is that the 
eigenfunctions do not concentrate at the singularity if the band index grows. 
A consequence is that the gaps are increasing and the transition matrix has 
the required properties. This result was announced in ||. In the second 
case for the countable set of resonant frequencies we prove the conjecture of 



18 1 concerning the location of the essential spectrum; a general argument 



based again on the off-diagonal decay of matrix elements allows to conclude 
that also in this case the spectrum is pure point. 

These results strongly suggest that in the models considered here, in fact 
for a < in the Ao language and for reasonable boundedness of the transi- 
tion matrix, pure pointness of the spectrum should not depend on number 
theoretical properties of the frequency! 

The fact that the spectrum of the Floquet operator is pure point- in the 
Hr picture- does not imply on general grounds that the energy expectation 
is bounded in time as the example in |14| shows. We prove here that applica- 
bility of the KAM method provides a uniform bound on the energy growth, 
so this applies here to the non- resonant case and holds true for a subclass of 
general time dependent Hamiltonians studied in p6[ |2^] which complements 
their results. 

The organization of the paper is as follows. In section ^| we fix notations, 
define the problem in detail and provide a regularization necessary for the 
methods in the following sections to work. Section |3| is devoted to the study 
of the non-resonant case; the result on pure pointness is Corollary FO, the 



bounds on energy growth are Corollary O and Theorem |3]5|. In section [| we 
determine the essential spectrum and prove pure pointness in the resonant 
case. 
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2 The problem and its Matrix Representa- 
tion 



We consider the class of potentials 

V{x) = P 5 '( x ~ nL ) + W ( x ) 

with W(x) = W(x + L) a differentiable multiplication operator on L 2 (M) and 
5' defined by the expressions ([!]) below. We refer to |§ and references therein 
for background material on this model. 

The limit (3 — > oo represents decoupling of the cells by a Neumann bound- 
ary condition. For (3 large the problem is a perturbation of the decoupled 
case, but only in the quadratic form sense. We shall show in this chapter 
that in spite of the singular character of the perturbation the problem can 
be represented by a matrix operator with polynomial off-diagonal decay. Let 
us fix 

Notations. L, F, (3 are positive numbers. We shall, however, mostly work 
with the parameters us := FL, g :— 1/(3 and employ the symbols T := J^- for 
the Bloch period; D := —id where d denotes a partial derivative; DM(OM) 
for the diagonal (off-diagonal) part of a matrix M; x f° r the binary code 
defined by x{True) := l,x(False) := 0; cte for a generic constant, indepen- 
dent of the parameters, which may change from line to line. H n denotes a 
Sobolev space of order n. We shall try to avoid to note the dependence of 
parameters of a quantity if we feel that this is possible while keeping clarity. 

The Stark- Wannier Hamiltonian is 

H S = D 2 X + Fx + W{x) 
defined on 

D(H S ) = G H 2 (R \ (ZL)), H s ^j G L 2 (R), 

ip'(nL+) - ip'(nL-) = 0, ip(nL+) - ip(nL-) = (3^'(nL)}. (1) 

The Hamiltonian for the driven ring is informally : 

H R (t) = V 2 + W(x) + (V5) on L^S 1 ), 

where V := D x — Ft; we use, however, the time dependent gauge transforma- 
tion exp (—iFtx) to transform the propagator into the periodic one of period 



4 



T generated by the Hamiltonian 

H{t) = H(t,u,g):=D 2 x + ^x + W(x) (2) 

D(H(t)) := {^Gif 2 ((0,L);e^y(L)-^(0) = 0, 
S (^(L)-#)) = ^'(0)}. 

Notice that because the domain is t dependent an argument for existence 
of the propagator is needed. This will be shown by mapping the problem to 



one whose propagator is known to exist, see remark |2.3| at the end of this 
section. 

We shall henceforth study the point spectrum of the Floquet Hamiltonian 

K = K(uj, g) = D\ er + H(t, u, g) (3) 

acting in L 2 ((0, T), dt; L 2 ((0, L), dx)) on the domain 

D(K) = {ip E H\(0, T), D(H(t, cu, g)), ip(T, x) = ^(0, x))}. 

An eigenvector of K with eigenvalue e will provide us with a Bloch- 
Floquet solution of the Schrodinger equation 

D t ip(t,x) + H(t)t/j(t,x) = 

which is of the form 

ip(t,x) = e~ iet (t){t,x) 

with <p periodic in t. 

A second reason to introduce K is the unitary equivalence of Hs and K, 
see i: 

U B H S = KU B 
where Ub is the Bloch transformation 



U B ■ L 2 {R) -> L 2 ((0,T) x (0,L),dt dx) 



T 

7« 
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The matrix representation M of K is constructed as follows: let {ip n (t)}neN 
be a periodic orthonormal eigenbasis of H(t): ip n {t + T) = ip n (t). 

{0i)ieZxN, <Pj(t,x) = <p jlij2 {t,x) := —=-if> ja (t,x) 

then is a basis of £ 2 ((0, T), dt; L 2 ((0, L), dx)). Using (., .) for the scalar prod- 
uct in x space we define 

M jk = «&,#0 fc )):=/ ^K<j> k ){t)dt (4) 

Jo 

T 

giCi-fciM (( fel£l , + E k2 (t))8 j2h2 + (^ 2 , AVU(*)) dt. 

In the rest of this section we shall study the properties of the eigenvalues 
E n (t,ui,g) of H and the coupling matrix (ip n , D t i() m ) . 

For ip G if 1 ((0, L)), <p G D(H(t)) we find by integration by parts: 

+ S (e^(i)-f(0))(e !U V(i) - 0(0)). 
So denoting the Neumann decoupled operator (g = 0) by 

H = D 2 x + yX + W(x) with ^'(0) = ip'(L) = 

we have the representation 

H = H + g\f(t,u))(f(t,u;)\ 

where 

f(t,u) = e-^5 L -5 . 

f is in if _1 ((0, L)) so we are in the framework of generalized rank-one pertur- 
bations; we shall use the results of ||27|| . H(t,u,g) is an analytic family with 
constant form domain -£^(0, L) for (t, u>, g) G S at xS aLJ xC for some a t , a u > 
where S a := {2; G C; |Imz| < a}. For the resolvent R(z) = (H — z)' 1 it holds: 

R(z) - R (z) = - 1+ 9 gG{z) \Mz)f)(Mz)f\ (5) 

with G(z) = G(z,t,u):= (f,Ro(z)f)- 

In the sequel we make statements for g small enough. This could be 
circumvented by the use of an adiabatic technique. We shall not do so as in 
section |3] the smallness of g will be essential anyhow. We obtain 
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Theorem 2.1 For g sufficiently small, u in a given interval [u_,u + ] C 
(0, oo), T = — , t G [0, T] £/ie operator H(t,u,g) as defined in equation (||) 
/ias simple discrete spectrum. For its eigenvalues E n = E n (t,u,g) it holds 
uniformly in t, u, g, n: 

(i) > cte > 0, 

(ii) < d u E n < cte < 1; 

furthermore there exists a basis {ipn} of eigenfunctions of H with ip n (t + T) = 
ip n (t), ip n G C°°([0, T] x [cj_,c<j + ] x [0,g m(M ;]) swc/i t/iat m i/ie C°° topology 
and uniformly in t,u,n: 

(Hi) (ip n , D t ip m ) = | ra fi^ 2 | /or fi/m, 

^ (if) n , D t if) n ) = 0(g/n). 

Proof. The behavior for large n has to be controlled. We compare 
H to the Neumann Laplacian -A^ on H 2 (0,L) with boundary conditions 
ip'(0) = ip'{L) = 0, whose eigenvalues are {im/L) 2 . This is done in two steps: 
first we compare H(t,u,0) -which is actually time independent- to — 
using regular perturbation theory; secondly the difference H(g) — H(g = 0) 
is treated using formula @. 

By a Wronskian argument the eigenvalues of H(t, u, 0) are simple and for 
n large enough it holds: 



The reason why the transition matrix decays and eigenvalues stay nearby 
upon switching on g is contained in the following auxiliary result. 

Lemma 2.2 For e L 2 ((0, L)) and the eigenprojections P n of H it holds in 
the C°° topology and uniformly in n: 



Proof. (Of the Lemma) By Riesz's formula we have with a circle T n of 
length |r„| centered at {im/L) 2 



E n (t,u,0) - (— ) 2 | < cte\\-x + W 



\(P n (t,u,g)<i>)(L)\ + \ (P n (t,u,g)<f))(0)\ < cte\\<j>\\ 



n 
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Denote be a, b indices which take the values and L. In order to prove 
the estimate on P n 0(a) = (5 a , P n (f>) by Krein's formula-equation ([5])- it is 
sufficient to show 

sup (\(8 a ,Ro(z)<f>)\ + \(8 a ,Ra(z)8 b )\) = 0( ' 



r, 



To do this we use regular perturbation theory and the fact that the "gaps" 



of the Neumann operator are growing. Denote d n := |4((n) 2 — {n — l) 2 ) 



and choose no such that sup w + W\\ < for n > n choose a suitable 
number M and and |r„| := Then it holds for the resolvent R N of the 



Neumann Laplacian 



M ' 



\\(%x + W)R N (z)\\ < — <l&nd \\R N (z)\\ = O(-) 
L n n 

uniformly for z G T n and 



Ro = R ^^ + ^ x+ w)R N r i 

= R N -R N (l + (^ rX + W)R N )- 1 (yX + W)R N . 

So the question is reduced to the explicit calculation of \\8 a R || and 
(8 a ,R N 8b)- R N (z) is given by its kernel R N (z)cj)(x) = J R N (x,y; z)4>(y) dy 

R N (x, y; z) : = 1 cos(yG(x A y)) cos{yfz{{y V x) - L)) 

y z sm zL 

with the notation x Ay (x V y)for the minimum (maximum) of x and y. One 
finds for example 

\\8 R N (z)\\ < || 1 cos^(y-L)\\=0(-) 

and similarly \\8 L R N (z)\\ = 0(±),\(5 a , R N (z)5 b )\ = 0(±) for a, b in {0,L} 
uniformly for z G T n . These estimates are preserved upon differentiation 
with respect to t, u, g. □ 



To continue with the proof of the theorem for the eigenvalues we show 
E n (t, lu, g) = (^) 2 + | + (W) + ^(1 - (-1)" cos^t) + 0(l/n) (6) 



in the C°° topology. (W) denotes the mean value (1/L) f Q W. Indeed by 
regular perturbation theory with the notation V := ^ + W + g\f)(f\ it is a 
corollary of the previous lemma that 

= — -trU (z - (^) 2 )R(z)VR N (z)VR N (z) dz) = 0(l/n) 
27TZ J Vn L 

With P^ = j \ cos ^x) (cos | (n > 1) the explicit term of the approxi- 
mation follows. From this we infer the assertions concerning the eigenvalues 
for n large enough. For the lowest finitely many n we employ continuity of 
d^En and compactness of to deduce (ii). 

The eigenfunctions are now constructed as 

P n {t,u,g)lP n {u) 

i/j n (t,uj,g) := — — 

\\P n {t,uj,g)ipv{uj)\\ 

for any time independent choice of eigenfunctions ip® of Hq. We differentiate 
HP = EP in the quadratic form sense to get P m d t P n = %%ffi ■ It follows 
from Lemma |2.2| for the off-diagonal part: 

(VwP^n) g(g) 

For the diagonal a calculation using Lemma TE?1 yields 



□ 



Remark 2.3 T/ie existence of the propagator U(t,s) of H(t) is a corol- 
lary of the preceding theorem: denote by J(t) the unitary between l 2 (N) 
and L 2 ((0,L)) which maps the n'th canonical base vector to ip n (t). Then 
J~ 1 (t)(D t + H(t))J(t) = D t + h(t) where the matrix operator h is defined by 

h is analytic with constant domain so its propagator u(t, s) exists. U is then 
given by 

U(t, s) = J(t)u(t, s)J~ 1 (s). 
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3 Stability for non-resonant frequencies 



In this section we shall employ the KAM algorithm to diagonalize the matrix 
M of the Floquet operator K. M is considered as a perturbation of its 
diagonal DM. For generic values of the frequency u the eigenvalues of DM 
form a dense subset of the real line ||16|1 . We shall show in Corollary |3l| 
that for a large set of "non-resonant" to the spectrum of K is pure point, 



in Corollary £374] and in Theorem |3.5| that the energy of the system stays 
bounded. 

In order to measure the decay of matrix elements consider the following 
Banach algebras (see |T§): let r, 5 > 0, C (0, oo), (x) := (1 + x 2 ) 1/2 , 

B(Q, r,5):={uj^ M{u) G B(l 2 (Z x N)); oo > ||M||n, r ,* := 
sup (l^M + l ^-^ l) 

The result of the KAM algorithm we need here is: 

Theorem 3.1 Let r G (0, oo) be large enough, Q = [u^,u + ] C (0, oo), 
M = M(u) = M*(uj) a family of matrix operators in l 2 (7h x N) such that 

||OM||n,o,T < oo; M jj = uji + e j2 (u) (j = (j u j 2 ) GZxN) with 

inf Ui „ > cte > 0, 



(e„-e m )(u)-(e n -e m )(i/) 



< 1. 



Then there is a 5 G (0, r) such that for 7 small enough and ||OM||n,o,r — 7 2 
there is a set of good frequencies C £1 with measure 

an<i a unitary family U^u) with || fToo ||r2 ; o,<s < °o suc/i t/iat 
U^MU' 1 ^) = M^H, OM oo = (cuGfioo). 

Remarks 3.2 

(%) In particular M(u) has a basis of eigenfunctions fj which decay polyno- 
mially: fj(k) = = 0(\k-j\- s ); 
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(ii) Actually 5 = r — cte so 5 can be chosen arbitrarily large if r is arbitrarily 
large. 



Outline of the Proof. The KAM method in its quantum guise first 
introduced by || is by now quite standard. We shall, however, give only a 
descriptive proof and refer to jy| and references therein for analytic details. 

The idea is to successively diminish the size of small off-diagonal elements 
by unitary transformations. We use the function x(Tr-ue) := l,x(False) : = 
0. Denote the matrices 

BMtj := M ijX (i = j), OM := M - DM, 

n 

D n M l - j := M ijX (\i -j\=n), B n M := £ D n M, 

and define recursively 

M. := BjM, U := I, 

W n;ij := X (i + 3), U n := e^U^, 

M n+1 := U n {B n+x M)U~ l = e w "M n e- w « + U n {T> n+l M)lJ- 1 . 

The idea of this is based on the identity 
ad Wn (DM n ) := [W n , DM n ] = -OM n 

which by the Lie Schwinger formula 

e w Me~ w = ^2 a< ^^ M ^ 

k=0 ^' 

leads to 

k 



n+1 = DM n + jr-Ty ad w n (OM n ) + U n B n+l MU-\ 
k=i ^ >' 



W n is to be estimated by OM n so the second term in OM„ +1 will be quadratic 
in ||OM n ||. It is in this estimate where one looses the resonant uo giving rise 
to small divisors. For each step one proves that for a > 1 and 7 n small 
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enough there is an open set fi n+1 C Q n with \Q n \ Q n+ i\ < cte fz~|j]^[j| such. 
that for r n+ i < r n it holds: 

cte 

\\W n \\n n+1>rn+ljS < — v2^+rl|OM||n n , r „,« I (7) 

fnV n i n+1) 

here |||e||| estimates the space part of DM„. The bad frequencies are con- 
trolled by a diophantine estimate 

^ n \^n+i= [J \uk + e m - e n \ < y(k + n - my 7 }. 

fc,m,n£N,n>m 

The growing gap property is then used to show that the contributions to 
the measure are summable. Estimating now ||ac?^M|| < cte ||W|| fc ||M|| it 
follows for r n+ i < r n with the shorthand || ■ || n := || • ||n„,r n ,5 



||OM n+1 || n+1 < (8) 
cte\\W4 n+1 e^ w »^\OM n \\ n + e 2 Xi l^""* 1 ||D n+1 M|| n , 

||DM n+1 -DM n || n+1 < 

cte\\W n \\ n+1 e^ Wn ^\\OM n \\ n + e 2 Z 3 ||D n+1 M|| n , 

\\V£\Ui < e* w ^\\U£.i\\n. 

The choice j n = 0(l/n^),r n = 0(l/n u ^ 1 ) for suitable n,v in estimate 
(||) then leads to a quadratic estimate for ||W n || n +i: 

||^n|U+i < c4W n \\l + ie cteWAn+1 +c 2 {nf 

where the constant c 2 is proportional to HOMH^o.r, (3 = 2// + v(2a + 1) — r. 
If ||OM||f},o,T is small enough and r large enough this implies that ||W„|| is 
summable, and that DM n and U n are convergent. □ 

As a consequence of this and the analysis in section |2| we obtain that the 
spectrum is pure point: 

Corollary 3.3 For the Floquet Hamiltonian defined in equation (Qj it holds: 

K(ou,g) has a basis of eigenvectors in L 2 ((0, T), dt; L 2 ((0, L), dx)) 

provided g is small enough and uj G Vt^ C [u;_,u; + ] the set constructed in 
Theorem \3. 1\ ) with measure \ [u^,uj + ] \ f^l = 0{^Jg). 
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Proof. Identifying L 2 ((0, L), dx) with / 2 (N) via the eigenbasis {ip n } of H 
constructed in Theorem |2.1| we find that K is unitarily equivalent to 

D t + E n (t)5 nm + (i/; n , D t ^ m ) 

on L 2 ((0, T), dt; / 2 (N)), which is of the form 

D t + h (uj,g) +gV(t,w,g) 

where V is C°° bounded and 0(1) in g, and ho(ui,g)ij := (Ej)5ij. Applying 
a version of the superadiabatic regularization as in [^T]],|T5], Thm.3.6] we 
get that K is unitarily equivalent to an operator of the same form whose 
fluctuating part has the property that (m 2 — n 2 ) T V mn is C°° bounded in I 2 
uniformly in the parameters for any r > 0. Furthermore, by Theorem |2| the 
diagonal elements still satisfy £n+ ^~ e " > cte > and < d u e n < 1 for g small 
enough. 

Finally going to the Fourier representation K turns out to be unitarily 
equivalent to a matrix M on Z 2 (Z x N) which has the properties required in 
Theorem □ 



By pure pointness of <j(K) every trajectory {t G K; U(t)ijj} generated by 
if is a precompact set. It follows that 

lim sup \\x(H(t) > r)U(t)if>\\ = 0. 

r ^°° tern 

This does, however, not imply that the energy expectation | (tp, H (t)ip) \ is 
bounded, as the example given in jnj shows. 

We shall show now that boundedness of the energy is in fact always 
ensured in cases where the KAM algorithm used above applies. 

Corollary 3.4 For g small enough there exists a set of frequencies C 
[u;_,u; + ] C (0, oo) with \[uj-,lu + ] \ Qoo\ = 0{- s /g) such that for uj G f2oo and 
for the propagator U of H(t, lu, g) it holds 

(1) U(t) = U p (t)e-^U-\Q) 

where G commutes and is relatively bounded with respect to H(0) and 
U p is T periodic and C°° as a bounded operator in L 2 ((0, L), dx); 

(2) \(U(t)iP,H(t)U(t)ilj}\ < cte 

for if) G Q(H(0)), uniformly for t G R. 
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Proof. Application of the KAM algorithm to the matrix M of K gave 
UooMU^ 1 = Mqo. We transform back to the space of (t,x) functions using 



the basis {cpj} whose space part is {ipj 2 } as constructed in Theorem |2TT 
Denote 

%o ■= {Uoo)jk\<Pj)(<Pk\- 

By construction Uoo is a Toeplitz matrix in the indices corresponding to the 
Fourier variable, i.e. (Uoo)jk is °f the form (U O0 )j 1 ^k 1 ,j 2 ,k 2 - Consequently 
is fibered, i.e.: 7^0 (t, a;) = (£)■?/> (t, x) for T^it) unitary, periodic and 
C°° bounded in L 2 ((0, L),dx); (M^)^ = uj x + ef 2 with ef 2 - (E j2 ) = 0(g) 
uniformly in n, uo. The reader may consult [15j for a more detailed discussion 
of this point. By Theorem [3J] this results in 

T^KT- 1 = £(MxOtf|&><&| = A + H^t) 

where is defined by (A + H 00 {t))'il) rn (t) = e^ipm(t). Denote 

E/ A (t) :=^|^ m (t))(^ m (0)| 

then the relation 

(A + H^ityUAit) = U A (t) ^e~P m (0) 
holds on D(H(0)) so with the definition G := ^e^P m (0) we obtain 

(U^T^HiT^UA) + ((t/ A 1 T 00 )A(T 00 1 f/ A )) = G. 

This formula implies the asserted form for the propagator with the defini- 
tion L/p(i) := T^U A {t) and the fact that U A is C°° bounded by Theorem |2.1|, 
furthermore it shows that U p preserves domains: U p (t)D(H(0)) C D(H(t)). 

For t/) in the form domain of H(0) it holds with tp := f/~ 1 (0)?/> 

(U(t)1>,H{t)U(t)i/>) = -(U(t)i>,D t U(t)il>) 
= (e- iG V, (G - (f/- 1 (Af/ P ))(t)e- tG V) 
= (^,GV> - (^,e^(C/ p - 1 (AC/ P ))(t)e- iG V) 
which is bounded uniformly in time as (?7_" 1 (A^»)) is periodic. □ 
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By the same method we complement now the results of [2fJ which were 



much extended in JJ^, [5|. These authors estimate the propagator by time 



dependent methods to discuss stability of the energy expectations. The spec- 
tral method used here is better suited to provide bounds valid on an infinite 
time scale. The differentiability properties on the potential could be relaxed. 
However, we do not make an effort, here, to do so. 

Theorem 3.5 Let T, g > 0, W a T periodic C°° function with values in the 
bounded operators on a Hilbert space. Consider H = ^2 n E n P n with growing 
gaps: n+ ^ a — - > cte > for an a > 0. 

Then it holds for the propagator U of H(t) := H Q + gW(t), ip G Q(H ): 

\{U{t)i),H{t)U{t)i)) \ < cte 

provided g is small enough and uj G fioo C C (0, oo), the set con- 

structed as in Theorem (\3. J\) with measure \[u_,u) + ] \ Q^l = 0(- s fg). 

Proof. By the KAM algorithm is applicable, so we can proceed as in 
the previous corollary and find U P ,G with U(t) = f/ p (t)e~* G *?7 J 7 1 (0). □ 



4 Stability for resonant frequencies 

In this section we shall show that for resonant frequencies to G Q(^) 2 the 
spectrum of K is still pure point. 

It is equivalent to show that the time T map U(T) of H has pure point 
spectrum. We go to the matrix representation of section ||. Let {ip n } be the 
basis found in Theorem pTT] , {e n } the standard basis of I 2 and J = J(t, cu, g) 
the unitary operator J : = J^o" \e n )(ip n \- We have 

J(t)(D t + H(t))J- 1 (t) = D t + h(t) 
h nm := E n 5 nm + (ip n , D t ip m ). 

An approximation of E n was worked out in equation ([]). Let g n (t) = 
g n (t + T) := J Q -^(— l) n+1 coscus ds, G the gauge transformation defined 
by G nm (t) = exp(ig n (t))5 nm . Now 

G(D t + h)G~ l - (A + ((^) 2 + | + (W) + f)6„ 
= 0(-)6 nm + e l ^- 9 ™\tp n ,D t ij m )- 
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by Theorem |2.1| the matrix function on the right hand side is a C°° function 
in the Hilbert Schmidt norm ||a||#s := (^2 nm \a n m\ 2 ) 1 ^ 2 and a forteriori in 
the compact operators on l 2 (N). We now make use of the argument of EnB 
and Veselic to conclude: 

Theorem 4.1 For u G Q(^) 2 ,g small enough, T = ^ it holds 

(i) a ess (U(T,co,g)) = (expH((^) 2 + f + (W) + f )T);n G N} ; 

L 2 ((0,L)) /ias a fraszs of eigenvectors of U(T,u, g) . 

Proof. Denote h := GhG' 1 + G(D t G~ l ), h : = ((^f) 2 + f + (W) + f ) X (n = 
m) and by [/, C/o their propagators. The spectrum 

a{U (T,u>,g)) = |ex P H((^) 2 + | + (W) + ^-)T),n G N 

is a discrete set. Furthermore 

J Uo\h-ho)U 

is compact for every s,t G R. By Theorem 5.2 of 01 ^"CO — Uq(T) is 
compact. So a ess (U(T)) = c ess (U (T)), which cannot contain continuous 
spectrum so a{U{T)) is pure point. G{T) = I, from the unitary equivalence 

U(T) = J-\T)U{p)J(T). 

we conclude that the spectrum of U(T) is pure point □ 



Acknowledgments 

Part of the results announced here was obtained during the visits of J. A. and 
P.D. to the Nuclear Physics Institute AS CR, Rez, and P.E. to the Universite 
de Toulon et du Var and CPT CNRS, Marseille; the authors express their 
gratitude to the hosts. The work was partially supported by the AS CR 
Grant No. 1048801 



16 



References 



[1] Ping Ao. Absence of localization in energy space of a Bloch electron 
driven by a constant electric force. Phys. Rev. B, 41:3998-4001, 1990. 

[2] J. Asch, P. Duclos, and P. Exner. Stark- Wannier hamiltonians with pure 
point spectrum. In M. Demuth and B.-W. , Schultze, editors, Partial 
Differential Operators and Mathematical Physics. Birkhauser, 1997. 

[3] J.E. Avron, P. Exner, and Y. Last. Periodic Schrodinger operators with 
large gaps and Wannier-Stark ladders. Phys. Rev. Lett., 72:896-899, 
1994. 

[4] J.E. Avron and J. Nemirovsky. Quasienergies, stark hamiltonians and 
energy growth for driven quantum rings. Phys. Rev. Lett., 68:2212-2215, 
1992. 

[5] J. Barbaroux and A. Joye. Expectation values of observables in time- 
dependent quantum mechanics. J. Stat. Phys., to appear, 1998. 

[6] J. Bellissard. Stability and instability in quantum mechanics. In S. Al- 
beverio and P. Blanchard, editors, Trends and Developments in the 
Eighties. Word Scientific, 1985. 

[7] F. Bentosela, R. Carmona, P. Duclos, B. Simon, B. Souillard, and 
R. Weder. Schrodinger operators with an electric field and random or 
deterministic potentials. Commun. Math. Phys., 88:387-397, 1983. 

[8] F. Bentosela, V. Grecchi, and F. Zironi. Approximate ladder of reso- 
nances in a semiinfinite crystal. J. Phys. C, 15:7119-7131, 1982. 

[9] A.M. Berezhkovski and A. A. Ovchinnikov. True width of electron levels 
in a crystal in a static electric field. Sov. Phys.: Solid State, 18:1908- 
1911, 1976. 

[10] L. Bunimovich, H.R. Jauslin, J.L. Lebowitz, A. Pellegrinotti, and 
P. Nielaba. Diffusive energy growth in classical and quantum driven 
oscillators. J.Stat.Phys., 62:793-817, 1991. 

[11] M. Combescure. The quantum stability problem for time-periodic per- 
turbations of the harmonic oscillator. Ann. Inst. H. Poincare, 47:62-82, 
451-454, 1987. 



17 



[12] S. De Bievre and G. Forni. Transport properties of kicked and quasi- 
periodic hamiltonians. Comm. Math. Phys., 1998, to appear. 

[13] C.R. de Oliveira. Some remarks concerning stability for nonstationary 
quantum systems. J. Stat. Phys., 78:1055-1066, 1995. 

[14] R. del Rio, S. Jitomirskaya, Y. Last, and B. Simon. Operators with 
singular continuous spectrum, iv. Hausdorff dimensions, rank one per- 
turbations, and localization. J. d'Analyse Math., 69:153-200, 1996. 

[15] P. Duclos and P. Stovicek. Floquet hamiltonians with pure point spec- 
trum. Commun. Math. Phys., 177:327-374, 1996. 

[16] P. Duclos, P. Stovicek, and M. Vittot. Pertubation of an Eigen-value 
from a dense Point Spectrum: a general Floquet hamiltonian. Preprint 
CPT, (P.3559), 1997. 

[17] V. Enss and K. Veselic. Bound states and propagating states for time- 
dependent hamiltonians. Ann. Inst. H.Poincare, sect. A, 39:159-191, 
1983. 

[18] P. Exner. The absence of the absolutely continuous spectrum for 5' 
Wannier-Stark ladders. J. Math. Phys., 36:4561-4570, 1995. 

[19] V. Grecchi, , and A. Sacchetti. Lifetime of the Wannier-Stark resonances 
and perturbation theory. Comm. Math. Phys., 185:359-378, 1993. 

[20] G. Hagedorn, M. Loss, and J. Slawny. Non stochasticity of time- 
dependent quadratic hamiltonians and spectra of canonical transforma- 
tions. J. Phys. A., 19:521-531, 1986. 

[21] J. Howland. Floquet operator with singular spectrum i,ii. 
Ann. Inst. H.Poincare, sect. A, 49:309-334, 1989. 

[22] A. Joye. Upper bounds for the energy expectation in time-dependent 
quantum mechanics. J. Stat. Phys., 85:575-606, 1996. 

[23] A. Kiselev. Some examples in one-dimensional geometric scattering on 
manifolds. J. Math. Anal. Appl, 212:263 ff., 1997. 



18 



[24] M. Maioli and A. Sacchetti. Absence of the absolutely continuous spec- 
trum for Stark-Bloch operators with strongly singular periodic poten- 
tials, erratum. J. Phys. A, 31:1115-1119, 1998. 

[25] E.E. Mendez and G. Bastard. Wannier Stark resonances. Physics Today, 
46(6):34-42, 1993. 

[26] G. Nenciu. Adiabatic theory: Stability of systems with increasing gaps. 
Ann. Inst.H.P., 67(4):411-424, 1997. 

[27] B. Simon. Spectral analysis of rank one perturbations and applications. 
In J. Feldman, R. Froese, and L. Rosen, editors, CRM Lecture Notes, 
volume 8. AMS, 1995. 



19 



